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Abstract 

A model for the supercooled liquid is considered by taking into account its solid like properties. 
We focus on how the long time dynamics is affected due to the coupling between the slowly 
decaying density fluctuations and the local displacement variables in the frozen liquid. Results from 
our model agree with the recent observation of Novikov and Sokolov [Nature (London) 431, 961 
(2004)] that the fragility index m of a glass forming material is linearly related to the corresponding 
ratio K/G of the bulk and the shear moduli. 







I. INTRODUCTION 



Understanding the basic mechanism for solidification of a supercooled liquid into an amor- 
phous structure constitutes an area of much current research interest in Condensed matter 
physics. As the liquid is increasingly supercooled its shear viscosity rj increases and the dy- 
namics slows down enormously. The glass transition temperature T g is characterized by the 
generic value of i] ~ 10 14 P with the corresponding relaxation time reaching the laboratory 
time scales. Experiments show that the nature of relaxation in different supercooled liquids 
approaching the glass transition is not universal. Liquids have been classified as fragile or 
strong depending on their dynamical behavior in the vicinity of T g . Strong liquids show a 
steady increase of viscosity with the lowering of temperature. In fragile systems the viscosity 
first increases slowly in the temperature range higher than T g and this trend is followed by 
a much sharper increase of r\ near T g . This classification of different glass forming materials 
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is facilitated in terms of the so called fragility parameter defined 

m = d\nr ] /d(T g /T)\ Tg . (1) 
More recently evidence of a dynamic crossover in the value of the viscosity induced by 



pressure change has also been reported [2j. An associated characteristic of the supercooled 
state is its solid like behavior [3J . This is manifested in the elastic response of the system to an 
applied stress. The elastic behavior persists for times shorter than the structural relaxation 
time. The frozen solid with amorphous structure has transverse sound modes in addition to 
the longitudinal modes which are present in the normal liquid state. Recently Sokolov and 
Novikov [4j have demonstrated that the fragility parameter m of a liquid is linearly related 
to the corresponding ratio K/G of its bulk and shear moduli, i.e., m — 17 = 29(K/G — 1). 
The dependence of the fragility and the vibrational properties of the liquid on the basic 
interaction potential was also tested recently with computer simulations {5]. This strengthens 
further the possibility of understanding the vibrational and relaxation properties of the 
frozen liquid from a common standpoint. 

In the supercooled state of the liquid, a tagged particle makes rattling motion and is tem- 
porarily trapped in the cage formed by its neighboring particles. The cage eventually breaks 
giving rise to continuous particle motion. The time for which the tagged particle is localized 
grows with supercooling. The time correlation function of the collective density fluctuations 
develops a plateau over intermediate times and eventually decays to zero. Ergodicity thus 



persists in the supercooled liquid over longest time scales. This behavior of the density cor- 
relation function is explained in the Mode coupling theory (MCT) [6] of supercooled liquids. 
The basic theoretical scheme used in this approach involves a formulation of the dynamics in 
terms of a few slow variables. In this paper we present using the mode coupling approach, a 
model which includes the solid like properties of a supercooled liquid. We study the effects 
on the dynamics due to the couplings of slowly decaying density fluctuations with the extra 



slow modes present in the amorphous so 
between m and K/G reported in Ref. 



id. Our results conform to the above linear relation 
4|. The paper is organized as follows. In the next 
section we present the model for the dynamics using an extended set of slow modes for the 
liquid state. In section III we present the results from the model and demonstrate how it 
explains the observed data for simple liquids. We end the paper with a discussion. 

II. MODEL STUDIED 

A. Nonlinear dynamics of Slow modes 

The dynamics of the liquid is formulated in terms of a small set of slow modes. The MCT 
is a basic step in this direction for studying the dynamics of a strongly interacting liquid (at 
high density) by taking into account the effects of coupling the slow modes in the system. 
The simplest version of MCT deals with the conserved densities of mass and momentum 
respectively denoted by p(x, t) and g(x, t). We assume the system to be a collection of N 
classical particles each of mass m. r a (t) and p* a (t) are respectively the position and the 
momentum of the a-th particle at time t. For the density p and the momentum density g 
we use the standard prescription 3], 



N 

p(x,t) = 5(x - f^it)) , 

a=l 

N 

&{2,t) = £*4(*W*-rS(*)) • (2) 

a=l 

For the solid like state in which the particles vibrate about their mean positions, the above 

fl 

set of slow modes is further extended to include the new collective variable |8j u(x, t). This 
is defined in terms of the displacements u a (t) of the a-th particle, (a = 1 to N) from their 
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respective mean position denoted by R° [9| such that R a (t) = R° + u a (t). We adopt the 
definition, 



N 



p(x, t)ui(x, t) =mJ2 - r «W) > (3) 

a=l 

The metastable positions of the atoms R° (a = 1 to N) in the glassy system remain 
unaltered for a long time. In a crystal with long range order these positions are independent 
of time. The equation of motion for u is obtained using standard procedure 0, [l(| in the 
form of a generalized Langevin equation, 

dui g gi 5F 

— + - ■ Vui = r *i^ + /» ( 4 ) 

at p p ouj 

where indices i,j refer to the different spatial coordinate axes. The Gaussian noise /j 
is related to the bare damping matrix T^-, through the fluctuation dissipation relation, 
< /i(x, t)/j(x', £') >= 2kBTTij5(x — x')S(t — t'), where ksT = (3~ l is the Boltzmann factor. 
The average thermal 

speed of a liquid particle of mass m is v — (3rh. F[p, g, u] is the effective Hamiltonian 
such that the probability of the equilibrium state is given by e~ F . For the isotropic solid F 
is obtained in the general form 



F= /- 
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(5) 

where 5p = p— p Q is the density fluctuation and p Q is the average mass density. The quantities 
A and B in ([5]) are related to the static structure factor ( correlation of density fluctuations 
at equal time ) for the amorphous solid. The symmetric strain tensor field is defined in 
terms of the gradients of the displacement field u(x), 2sij = (V%Uj + Vj«j) — V 'iU m V jU m 
( such that Sij = Sji). The trace and the traceless parts of Sy, respectively defined as 
st = J2i s a an d Sij = Sij — 5ijST/3, appear in the expression (0) for free energy of the 
isotropic solid. 

The equation of motion for p is the continuity equation, 

! + V. g = 0. (6) 

For the momentum density g the dynamics is given by the generalized Navier-Stokes equa- 
tion, 

| + Ev^ = ^, (7) 



where 0j's are the gaussian noises related to the bare or short time transport properties of 
the liquid jl2l ]. The symmetric stress tensor cr^- is obtained as a sum of the reversible part 
afj and the irreversible (dissipative) part a^, such that {9} 

4= ?f +Piti - 2G s il + 2s iA Ks T + BS A+ ^s im (8) 

where K = K — 2G/3. The quantity P in (|SJ) is identified with pressure in conventional 
hydrodynamics and is obtained as local functional of the hydrodynamic fluctuations, 

X 2 2 

P = (A Po - B)5p + (B Po - K)s T + A^- Gs, mSnrf . (9) 

The dissipative part of the stress tensor <r^ is expressed in terms of the bare viscosities. 

= -VoFiVj + Vj«i - 3<%(V.v)] - Co%(V.v) (10) 
where rjo and Co are the bare shear and bare bulk viscosities respectively and v = g/p. 



B. Renormalization of transport coefficients 

The time correlation functions of the slow modes are obtained by averaging over the 



noises in the nonlinear equations of motion using standard field theoretic methods 13] . The 
correlation functions between two slow modes ipa and ipp are defined as, 



G aP {\2) =< MVM 2 ) > (ii) 

The dispersion relations for the various hydrodynamic modes in the system are obtained 
from the pole structures of the correlation functions. The corrections to the correlation 
functions due to the nonlinear it ies in the equations of motion for the slow variables are 
expressed in terms of the self-energy E defined through the Dyson equation 

G- 1 (q,^) = [G°(q,c)]- 1 -S(q,^). (12) 

G°(q, uj) refers to the matrix of the correlation functions corresponding to the linear dynam- 
ics of the fluctuations. Assuming the amorphous solid state to be isotropic the correlation 
function can be separated in terms of a longitudinal and a transverse part as 

G<Hto u ) = MjG^iq, u) + (5 i:j - q4j)G T a ^ u) (13) 
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Similarly the self-energies are split into their respective longitudinal and transverse parts 
denoted by and Sy. Using the Dyson equation ffl2|) the renormalized viscosity of the 
liquid taking into account the nonlinearities in the dynamics is obtained in the form 

i L 9 M^) ( 14 ) 

where 



r 
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s se (q^) = -q 7, 
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(15) 

The functional forms of the self energies are computed with the diagrammatic techniques 
of the MSR field theory JjJ. Of particular interest in the context of glassy dynamics is 
the density auto correlation function <f)(q,t). The Laplace transform of 4>(q,t) (which is 
normalized with respect to its equal time value) is obtained as a two step continued fraction 
in terms of the memory function T(q, z) 12], 

1 



<f>{q,z) 



(16) 



Qq is the microscopic frequency of the liquid state[l5J]. For the case of a normal liquid with 
only the standard set of conserved densities as a slow variable, the dominant nonlinearity in 
the momentum equation comes from the coupling of the density fluctuations in the Pressure 
functional in (jSJ). This obtains 16] the standard mode coupling model in which memory 
function includes only products of density correlation functions 17L Il8j|. The present formu- 
lation involving an extended set of hydrodynamic variables which include u ( in addition to 
the conserved densities) obtains the longitudinal as well as the transverse sound modes in 
the amorphous solid. Their speeds are respectively obtained as, c 2 L = Mj p Q + Ap — 2B and 



4 = G/p , where M = K + 4G/3 is the longitudinal modulus. By including the coupling 
of density fluctuations Sp with the displacement field u ( see eqn. (JH])), a new contribution 
to the memory function T(q,t) is obtained. At the one loop level this is obtained as, 



r(q,^) = r + r( 1 )(q,^) + r( 2 )(q,^) 

rW(q,w) and r^ 2 ^(q, lo) are given by the expressions 



(17) 



(18) 
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and 



r<2 ' < q ' ^ = kjf I W / ~ k ' " " < 19 ) 

The role of convective nonlinearities is assumed to be absorbed in T which is the bare 
part or short time part of the transport coefficients. The contribution is obtained from 
the first diagram of fig. 1 resulting from the coupling of the displacement field u with the 
density fluctuation Sp. is the contribution from the second diagram in fig. 1 coming 
from the coupling of density fluctuations. In evaluating these diagrammatic contributions we 
make the approximation that for time scales ( short compared to the structural relaxations) 
over which the supercooled liquid displays elastic behavior, the longitudinal and transverse 
correlations of the displacement field are frozen ( constant in time ), i.e., ~ 5 (us). 

Therefore as a result of the solid like behavior over intermediate time scales, G^ u and G^ u 
are then obtained as 

GL(Kt) = j§, G T uu (k,t) = ^ (20) 

where the k~ 2 dependence of the correlation function arises from the static structure factor 
for the displacement fields. Substituting these approximate forms for the displacement cor- 
relations and evaluating the integrals in ([18]) and (JT9J) , we obtain for the density correlation 
function and the memory function the simple form of coupled nonlinear integral equations, 



= .-»;/[, + r W1 < 21 » 
r(t) = Cl 0OO + c 2 2 (t) (22) 

where in eqns. (1211) and (|22|) we have dropped the wave vector dependence of and V 
for simplicity. The density correlation function is approximated in the form G pp (q,t) = 
Xp P (g)0(t) with Xpp being the equal time correlation of the density correlation function 
determined by the thermodynamic properties like temperature and density. The constants c\ 
and C2 in (122j) are determined from an evaluation of the vertex functions in this approximation 
( of wave vector independence) as, 

ci = 2 + /(*)], c 2 = A ° 2 , (23) 
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where A CT = A (l — 2c)/(2 — 2a), and /(cr) = 2/5(1 — 2a) with a being the Poisson's 
ratio a = (3K — 2G)/[2(3K + G)]. We have used in the expressions fl23|) the definitions 
A G = B 2 /(AG) and A Q = (A 3 /67r 2 72 )(t>/cL) 2 in terms of parameters dependent on the 
thermodynamic state of the system. The equilibrium number density of particles in the 
fluid is n (p D = mn ). A is the upper cutoff of wave vector representing the shortest length 
up to which the fluctuations are considered. In considering the mode coupling expression 
for the viscosity we have ignored here the presence of very slow vacancy diffusion mode and 
its coupling to density fluctuations. Finally, it is useful to note here that we are considering 
the simple form of the model in which all processes giving rise to ergodic behavior in the 
asymptotic dynamics have been ignored. 



III. RESULTS 



The central focus in the present analysis is the time dependence of the density autocorre- 
lation function <p{t). From the coupled set of eqns. ( fl6l) (inverse Laplace transformed in the 
time space) and (1221) . we obtain a nonlinear integro-differential equation for the dynamics 
of cf>(t), 

4>{t) + 000 + 0(t) + f dsT[<p{t - s)]4>{s) = (24) 



where the dots refer to derivative with respect to time. This equation is solved numerically 
to obtain the behavior of the density correlation function. It is clear that the dynamics of 
the density correlation function <f)(t) is driven here by the memory function T(t) which is 
expressed in terms of itself. This constitutes a nonlinear feed back mechanism [itI . fis| 
causing a dynamic transition of the liquid to a nonergodic state in which the long time limit 
of the density correlation function cf)(t — > oo) = /. The quantity / is termed as the non 
ergodicity parameter (NEP). In the plane of C\ and c 2 the dynamic transition occurs along 
the line ci = 2 v /c 2 ~ - c 2 [19|]. 

In fig. 2 we display the phase diagram with c\ and c 2 showing the ideal glass transition 
line. Though this dynamic transition is finally cutoff due to ergodicity restoring mechanisms, 
it marks a cross over in the dynamical behavior of supercooled liquid. The liquid state is 
characterized by the density correlation function following initial power law relaxations, 
followed by a final stretched exponential decay ~ exp [-{t/r) } [19|, |20|. In the glassy 



state the density correlation function decays to a nonzero value equal to the non-ergodicity 
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parameter. In fig. 2 we have also shown the values of c\ and C2 corresponding to the 
nonergodic states in which the dynamics is studied in this paper (to be explained below). 
For presenting our results in the following, we express time in units of V j {p Cl) where 
r o = Co + 4r/ /3 is the bare longitudinal viscosity of the liquid related to its short time 
dynamics. We keep the A c appearing in the expression (I23p for the coupling constants of the 
memory function fixed at a constant value (=.4) throughout the calculation. This numerical 
value is reached by treating A as an adjustable parameter here for comparison of the results 
of the present model with experimental data. This essentially implies that the cutoff A of 
wave vector integration is being treated as an adjustable parameter in the coarse grained 
model we present here. 

A key quantity of interest in the present analysis is the fragility parameter m which 
by definition is related to the final relaxation behavior of the supercooled liquid near T g . 
But over this temperature range, the MOT approach has uot beeu very successful in 
explaining the relaxation behavior. A useful observation [4] in this respect is that the slope 
of the viscosity vs. inverse-temperature plot in the high temperature range can be linked 
with the fragility. This is justified as follows : For very fragile systems (large m) the slope of 
the f] vs. T g /T curve at temperatures near T g is large and hence it must be correspondingly 
small at the other end of the Angell plot, i.e., for temperatures much higher than T g . This 
is because the curves for different m meet at both ends on the T g /T axis. By examining 
experimental data, it was pointed out [4] that the slope of the r\ vs. T g /T curve on the 
high temperature side is inversely related to the fragility m. Since MCT is a valid theory 
for the slow dynamics well above T g , we have investigated the relaxation behavior which 
follows from the present model in this high temperature range. We focus on the growth 
of the relaxation time r instead of that of the viscosity rj. In our model, increasing the 
parameter A brings the system closer to the ideal transition ( since it results in an increase 
of t ) and hence this parameter is treated like the inverse of temperature, A D oc T g /T. The 
dependence of the relaxation time r on A G at various (fixed) Kj G values is displayed in fig. 
3. For large K/G the observed variation is similar to the r\ vs. T g /T curve of a fragile liquid. 
Although according to MCT model the power law behavior describes the relaxation time 
data in this range better, we fit r, following Ref. [4J, to an Arrhenius form with activation 
energy k, i.e., r ~ exp(fi;A ). By the property of the Angell plot referred to above, we 
assume fh — 1/k to be proportional to the fragility m. A D is proportional to T g /T with 
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the proportionality constant being independent of K/G. Thus we assume that constants A 
and B in the expression (jSJ) for the free energy are only functions of temperature and the 
two elastic constants K and G are to be treated as independent entities. Fig. 4 displays 
m against the variation of K/G showing an approximately linear behavior. Therefore the 
present model agrees qualitatively with the linear correlation predicted by Novikov and 
Sokolov. It should be noted that the inverse nature of the relations between fragility and 
the slope of the viscosity vs. inverse-temperature plot respectively at the high and low 
temperature sides is a property of the Angell plot itself. However the reciprocal relation 
i.e. k oc I/771 as used above can only be justified in a qualitative manner. In order to 
demonstrate the sensitivity of the results on the values of the parameters A Q and A Q we 
display in fig 4. the m vs K/G curves for different values of Ao ranging from Ao = 0.3 to 
Ao = 0.8. The nature of the curves are qualitatively similar to what the best fit value for A D 
obtains. 

Next we consider the dynamics on the other side of the transition when the density 
correlation function freezes to a constant until ergodicity restoring processes take over. In 
the glassy state, the NEP / is estimated from the plateau which <p(t) reaches over long time. 
Fig. 5 displays the variation of / with the ratio K/G. The nonergodic behavior and the 
elastic response of the supercooled liquid as seen here are valid over initial time scales for 
which the system is frozen and is solid like. However estimating the fragility parameter m 
directly involves computing the temperature dependence of viscosity or relaxation time near 
T g . As already pointed out, the long time dynamics in the deeply supercooled state close to 
T g is still beyond the scope of MCT. Our calculation of the NEP at low temperatures here 
only refer to intermediate time scales up to which the present description in terms of the 
supercooled liquid is valid. Therefore in order to link fragility m with the elastic properties 
of the supercooled liquid we make use of the dependence of the fragility m on the NEP / 
as inferred from analysis of experimental data 211 ] . In fig. 6 the data points and best fit 
curve for the experimental results linking NEP / with the fragility m, i.e., m = 167 — 176/, 
is shown. The NEP values noted here are taken at the glass transition temperature T g 



21j. To test this empirical relation further we compare the a a = (1 — /)// vs. m in the 



inset of fig. 6 with another set of experimental data 22j from X-ray Brillouin scattering. 
Similar qualitative behavior is seen validating the link between the NEP with the fragility. 
The dependence of m on the ratio K/G of elastic constants is then obtained through their 
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common dependence on the NEP /. The variation of m with Kj G is shown in fig. 7. We also 
display the experimental data (from which the linear relation of Ref. 2jj was proposed) in 
the same figure, showing reasonable agreement with the predictions of the present model. In 
reaching this agreement between the theoretical model and experimental data we adjust the 
temperature dependent parameter A D to the value 1.7. As argued above A Q is proportional 
to Tg/T and hence if we assume near T g , A D = c (T g /T), then c Q = 1.7 obtains the best fit 
of the present model's predictions with results of Ref. 4j. To demonstrate the sensitivity of 
the results on the values of the adjusted parameters we display in fig 7. ( as in fig 4 for rh 
) the fragility m vs K/G curve for different values of Ao ranging from Ao = 0.3 to Ao = 0.8. 
A value is kept same for all the curves i.e. A = 1.7. Usually fragility index measurements 
refer to the slope of the viscosity vs. inverse temperature plot at the calorimetric glass 
transition temperature T g . Hence the parameter A D has been kept fixed here. The roughly 
linear behavior of the m vs. K/G curve is retained over the whole range of the value of 
the parameter Ao- The variations of / and m with K/G as displayed in the fig. 5 and 7 
respectively are obtained at this same fixed value for A G . With the variation of K/G, the 
corresponding values of c± and C2 are changed as shown by filled circles in the fig. 2. With 
increasing value of K/G we move to regions of small values of the couplings C\ and 0%. 

Finally we make two predictions using the proposed model linking the elastic behavior 
with relaxation properties in the glassy systems. According to the present model the density 
correlation function <p(t) has a complex relaxation behavior. On the ergodic side of the 
transition line, the initial power law relaxations change over to a stretched exponential 
decay with a stretching exponent (3. As a further illustration of the present model, we 
display in fig. 8 the variation of the stretching exponent (3 with the Poisson's ratio o = 
(3K — 2G) /[2(3K + G)]. For the results shown here we keep A D fixed at a typical liquid state 
value A = 0.885. It should be noted that following the discussion above the j3 presented 
here corresponds to temperatures much above T g . On the other hand for computing the 
fragility m vs. Poisson's ratio a in terms of the non ergodicity parameter /, as described 
above in Fig. 5-7, we work on the nonergodic side of the transition line , i. e. in the glassy 
phase where the freezing of the density correlation function is clearly visible (A G = 1.7). 
We do not have in the present mode coupling model a way of computing the (3 in this low 
temperature regime which will require taking into account the role of ergodicity restoring 
mechanisms in the long time dynamics. Hence the predicted behavior of the stretching 
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exponent should be tested at higher temperatures away from T g . 

Another suitable quantity is the power law exponent of relaxation near T g . On the non- 
ergodic side of transition line, the system is dynamically frozen initially until the ergodicity 
restoring processes take over. For initial time scales ( beyond the microscopic times) the 
density correlation function develops a plateau / with a power law form <f>(t) = f + Ct~ a 
decay. From the present model we are able to predict the dependence of the power law 
exponent a on the K/G ratio. Fig. 9 displays the dependence of the exponent a on K/G. 
The inset shows the nature of the power law behavior at a particular value of K/G = 2.6. 
The curve is obtained at Ao = 1.7. According to the predictions of the present model the 
power law exponent grows sharply for less fragile systems indicating a sharper freezing of 
the system over intermediate time scales for strong liquids. 



IV. DISCUSSION 



The present model for the amorphous system is constructed in terms of an extended set 
of slow modes. The approximations which go into the present formulation of the dynamics 
for the solid like state are well known and standard. In order to define the slow variable 
which accounts for the solid like nature of the amorphous state, we require reference to a 
rigid lattice. On the other hand the ergodicity restoring process in the system invalidates 
the existence of any such rigid structure over longest time scales. Therefore the elastic 
description of the liquid is only valid up to the time scale of structural relaxation. A self- 
consistent treatment with dynamic connections between the elastic and viscous behavior of 
the system will possibly provide a way to make a unified description including both of the 
above two aspects of the glassy state. Furthermore, the atoms in the solid vibrate around 
the random lattice sites over intermediate time scales. Generally strong glassy systems 



231 ] . In the latter 



have a higher network forming tendency than the highly fragile systems 
structural degradation occurs more easily. Thus for systems with very high m values the 
underlying assumptions of the present model hold for shorter time scales. 

The present model which considers the relation of the fragility of the liquid to its elastic 
properties brings out some basic characteristics of the dynamics. In fig. 7 we se that the 
increase in the fragility m corresponds to higher values of the ratio K/G while from fig. 2 
it follows that coupling c\ and C2 between the slow modes is weaker for larger values of this 
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ratio. Therefore the trend of the experimental data indicates that in a more fragile system 
there is a weaker coupling of the hydrodynamic modes. This is particularly relevant since 
the coupling between the slow modes constitutes ( through the formulation of the mode 
coupling theories ) the basic mechanism for slow dynamics in supercooled liquids. Also, 
since it follows from the present theory that larger values of K/G give rise to lower values 
of the nonergodicity parameter, the strength of the freezing or jamming of the system is 
weaker in more fragile systems. 

It is important to note that the present model is formulated for the case of simple liquids 
and accordingly it is expected to apply to only specific types of systems. The relation between 
m and Kj G stated by Sokolov and Novikov is based on the data analysis of non metallic and 
non -polymeric systems [24]. Efforts to apply this scaling to a wider group including metallic 



glasses or multi component systems fail 



25| for obvious reasons. Polymeric materials also 



clearly deviate from this linear relation. Monomers follow similar correlations, but increase 

in the chain length in some polymers leads to the deviations [26j. 

The model similar to the one used here having the memory kernel with both linear and 

□ 

quadratic terms has been initially proposed in Ref. [20(. However the present form of 
the r(t) ( having direct relation to the elastic constants) is only justified by including the 
displacement variable for the amorphous solid in the theoretical description. In this 
regard it is useful to note that the past literature on MCT can be broadly divided into two 
different groups. In one approach the dynamics of a supercooled liquid is studied starting 
from the basic equations which apply to the liquid state or crystalline state of matter. This 
is the scheme we have followed in our analysis here. The set of nonlinear equations we have 
considered for constructing the mode coupling equations in the present context are a) the 
conservation laws for the mass and momentum and b) the dynamical equation for the extra 
slow modes for the solid like state. The resulting model linking the non ergodic behavior 
to the elastic properties then follows through a careful consideration of the nonlinearities in 
the dynamics. The present model thus extends the standard MCT in a physically relevant 
manner to include the solid like nature of the supercooled state. On the other hand the MCT 
literature also have schematic models in which ad-hoc forms of memory functions are simply 
assumed in order to consider different types of possible relaxation scenarios which follow from 
the nonlinear feed back mechanism. It is only in this respect the present model is technically 
similar to the schematic 0i2 MCT studied by Gotze. Finally, we have ignored here the wave 
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vector dependence of the model to focus on the basic feed back mechanism. The parameter 
A G has been treated as a temperature dependent parameter ( ~ T g /T ) in the model. The 
temperature and density dependence of the results should follow in a natural way when 
such extensions of the model are considered. This calculation will require implementing the 
proper structure factor of the liquid in the formulation and will be considered elsewhere. 
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FIG. 1: Relevent one loop diagrams which contribute to the effective viscosity as non- linear 
corrections. 
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FIG. 2: ci vs. C2 and the solid line displays the dynamic transition line. The circles (joined with 
solid line to show continuity) shown in the nonergodic phase correspond to the states for which 
the results of fig. 5 and 7 are obtained (Ao = 1.7). These circles correspond to different values 
of the ratio K/G of bulk and shear moduli with the adjacent number to a circle indicating the 
corresponding value of K/G. 
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FIG. 3: Log of relaxation time r in units of To/(p cj j ), vs Ao (see text). The constant K/G value 
corresponding to each curve is also displayed adjacent to it in the figure. The dashed and dotted 
lines are, respectively, the power law and exponential fits to the corresponding curve. 
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FIG. 5: Non ergodicity parameter /vs the ratio K/G of bulk and shear moduli from the present 
model (A = 1.7). 




FIG. 6: Experimental data (points) and best fit (line) to the experimental data Ref. 
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fragility index m with /. The inset shows a Q = (1 — /)// vs fragility index m{ solid line) (using 



the experimental best fit), X-ray Brillouin scattering data ( solid triangles) 
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FIG. 7: Fragility index m vs ratio K/G of bulk and shear moduli from the present model at 
different Aq values (solid lines) (Aq = 1-7), experimental data of Ref. [4] (points). 
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FIG. 8: Stretched exponential exponent (3 vs Poisson's ratio a (Aq = 0.885). 
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FIG. 9: Power law exponent a vs the ratio K/G of bulk and shear moduli (Ao = 1.7). Inset shows 
normalized density correlation function (f)(t) vs t in units of ro/(/9 c| / )(see text). The dashed curve 
is the power law fit with the exponent value shown as a dot on the main figure. 



19 



